Conceptual Model

Model Implementation

Load Libraries

library(readr)
library(janitor)

Load Model and Data

# Load Model
source("almond_yield_anomoly.R")

# Load Data
data <- read_delim("clim.txt", delim = " ", col_names = TRUE, quote = "") %>%
  # Clean column names
  clean_names()

Run Model

# Run Model
almond_yield_anomaly_from_daily(data)
## # A tibble: 22 × 3
##       wy   Tn2     P1
##    <dbl> <dbl>  <dbl>
##  1  1989  8.64   2.80
##  2  1990  8.68  55.8 
##  3  1991 10.4  135.  
##  4  1992 11.9   69.6 
##  5  1993 10.9   77.9 
##  6  1994  9.62  34.8 
##  7  1995 12.6  677.  
##  8  1996 12.0   40.3 
##  9  1997  9.99 285.  
## 10  1998 11.5   89.8 
## # ℹ 12 more rows
## $yield_anomalies
##         1989         1990         1991         1992         1993         1994 
##   -0.3552237    9.2906757   68.9130633   15.4280698   20.2083803    2.4820009 
##         1995         1996         1997         1998         1999         2000 
## 1919.9811511    3.5818399  329.6938750   27.8636956   -0.1436364    9.5999883 
##         2001         2002         2003         2004         2005         2006 
##  159.5119587    0.2450914   -0.2585997   -0.2367722  656.3724121   18.6324135 
##         2007         2008         2009         2010 
##   20.2007396  576.2821943    0.7367438  153.7655092 
## 
## $mean_yield_anomaly
## [1] 181.4453
## 
## $max_yield_anomaly
## [1] 1919.981
## 
## $min_yield_anomaly
## [1] -0.3552237
# Save results
almond <- almond_yield_anomaly_from_daily(data)
## # A tibble: 22 × 3
##       wy   Tn2     P1
##    <dbl> <dbl>  <dbl>
##  1  1989  8.64   2.80
##  2  1990  8.68  55.8 
##  3  1991 10.4  135.  
##  4  1992 11.9   69.6 
##  5  1993 10.9   77.9 
##  6  1994  9.62  34.8 
##  7  1995 12.6  677.  
##  8  1996 12.0   40.3 
##  9  1997  9.99 285.  
## 10  1998 11.5   89.8 
## # ℹ 12 more rows
# View results (in tons/acre)
almond$min_yield_anomaly
## [1] -0.3552237
almond$mean_yield_anomaly
## [1] 181.4453
almond$max_yield_anomaly
## [1] 1919.981